Nonequilibrium quasiparticles and 2e periodicity in single-Cooper-pair transistors* 
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We have fabricated single-Cooper-pair transistors in which the spatial profile of the supercon- 
ducting gap energy was controlled by oxygen doping. The profile dramatically affects the switching 
current vs. gate voltage curve of the transistor, changing its period from le to 2e. A model based on 
nonequilibrium quasiparticles in the leads explains our results, including the surprising observation 
that even devices with a clean 2e period are "poisoned" by small numbers of these quasiparticles. 



Coherent superpositions of charge states can be pre- 
pared and manipulated in circuits made of ultrasmall su- 
perconducting junctions 0, 0. This phenomenon may 
enable the construction of solid-state qubits based on 
charge states, as well as a quantum current standard 
whose speed is not limited by the stochastic nature of in- 
coherent tunneling. In both cases, tunneling of unpaired 
quasiparticles (QPs) causes decoherence and may limit 
operation to impractically short timescales. Reducing 
this "QP poisoning" requires a detailed understanding 
at a fundamental level. 

The single-Cooper-pair transistor (SCPT), shown in 
Figure GJa), consists of a micrometer-sized island that 
has a capacitive gate electrode and is probed by two 
Josephson junctions with areas ~(100 nm) 2 . The island 
charging energy is modulated by the gate according to 
E c( n g) = E co(n - n g ) 2 , where E C o = e 2 /2C s , e is 
the electron charge, is the total island capacitance, 
rig = CgVg/e is the normalized gate polarization, and 
n is the integer number of excess charges on the island. 
This Coulomb energy suppresses fluctuations in the is- 
land charge and causes the Josephson energy Ej to vary 
with rig through charge-phase duality |2|. This effect is 
strongest when Eco > Ej ^> fcT, where k is the Boltz- 
mann constant and T is the temperature. In a current- 
biased configuration, the modulation of Ej manifests it- 
self in the current I sw at which the SCPT switches from 
the "supercurrent branch" near zero voltage to the "volt- 
age state" 4]. I swing) is maximized when charge states 
differing by one Cooper pair are degenerate, yielding a 
2e periodic modulation in n g . Without QP tunneling, 
each 2e interval of I sw (n g ) has a single peak, decreas- 
ing monotonically to a valley on either side. This shape, 
which we call "clean" 2e modulation, is illustrated in Fig. 
da) for parity labeled "even" or "odd" depending on 
n. Intermittent tunneling of QPs during a measurement 
causes random changes in parity and results in a more 
complicated modulation containing secondary peaks and 
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FIG. 1: (a) Schematic representation of Itw en • Single- 
Cooper-pair transistor circuit shown at right. (b)-(c) 
Grayscale plots of I sw i n g) f° r co-deposited type L and type H 
SCPTs at T = 30 mK, measured with J = 1 /xA/s. The cor- 
responding gap profiles are shown on the right. 



valleys. Thus clean 2e modulation has often been viewed 
as an indication that an SCPT is free of QP poisoning. 

Early studies of SCPTs (many unpublished) showed 
a le period, indicating strong QP poisoning, even when 
thermal QPs were suppressed. Nonequilibrium QPs in 
the leads were thought to be responsible, leading one 
group to put normal metal leads close to their junctions 
to filter QPs |3(. This yielded a clean 2e period in almost 
all devices and appeared to resolve the issue. However, 
other experiments revealed a variety of unexplained re- 
sults, including a 2e period without QP filters [1,0, 2e 
in some devices and le in others even though all devices 
shared the same design, fabrication process, and mea- 
surement setup > an d finally poisoning that increased 
as T was lowered in some samples 9]. In this sense, QP 
poisoning in SCPTs has remained a puzzle, which moti- 
vated us to create devices that were controllably either 
le or 2e and understand in detail how QPs behave in 
both types. 

SCPT fabrication typically involves a double-angle Al 
deposition, with a junction oxidation step between lay- 
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ers, that creates the leads in one deposition and the is- 
land in the other. Since the process is completed without 
breaking vacuum, the Al films are usually assumed to be 
identical, but minor variations in vacuum conditions may 
cause differences in film disorder that can affect the su- 
perconducting gap energy A. Thus an SCPT may have 
a lead gap A^ that is different from the island gap A^. 
We have exaggerated this effect by oxygen-doping the 
first deposition allowing us to co- fabricate SCPTs 
in which 5 A = Ai — Ai is both positive (type L) and 
negative (type H) (see Fig. Elb) and (c)). Starting from 
a base pressure ^3xl0 -7 mbar, we deposited the first 
Al layer (20 nm thick, 0.4 nm/s) while flowing O2 gas 
to raise the pressure to ^5xl0 -6 mbar. We then oxi- 
dized at room temperature in 130 mbar of O2 for 5 min- 
utes. After rotating the substrate to the second angle and 
pumping back to base pressure, we deposited the second 
layer (30 nm thick, 0.4 nm/s) with no O2 added. Mea- 
surements of individual films showed the oxygen-doped 
deposition always had a higher critical temperature T c 
and thus a larger A. We note that our devices have no 
QP filters near the junctions, however the Al leads con- 
nect to Au/Ti pads c± 10 jam away. 

All measurements were performed in a dilution refrig- 
erator with a base temperature of 25 mK using coaxial 
leads with microwave filtering at both 4 K and the mixing 
chamber. The devices were mounted inside an rf-tight 
copper box and measured using a two-probe, current- 
biased configuration. The bias current was ramped using 
a linear voltage ramp applied across a 10 resistor 
(at 4 K) in series with the SCPT. In this configuration 
the current-voltage characteristic was hysteretic with a 
distinct super current branch at zero bias. We measured 
I 3W , the current at which the SCPT switched to the volt- 
age state, by cycling through the hysteresis loop 10 3 -10 4 
times at each value of gate voltage. After each switch, 
we returned to the supercurrent branch and allowed the 
system to equilibrate for at least 1 ms before starting the 
next ramp (longer times did not change our results). In 
this manner, we acquired a histogram of I sw at each value 
of n g , allowing us to study the distribution of switching 
currents rather than only the mean as in most previous 
studies of SCPTs. 

Here we focus on two co-deposited devices which are 
representative of the L and H gap profiles shown in Fig.Q 
and which are otherwise nearly identical. We deter- 
mined the gap energies of the first and second depo- 
sitions from the measured critical temperatures, yield- 
ing Ai = 246 fieV (T cl = 1.63 K) and A 2 = 205 fieV 
(T C 2 = 1.36 K). We determined the charging energies 
using the asymptotic current- voltage characteristics [ll| , 
giving Eco — 115 jaeV for both devices, while the to- 
tal normal state resistances, Rn,l — 19 k£l and Rn,h = 
18 kO, gave us Ambegaokar-Baratoff values for the 
Josephson energies per junction of Ej^ = 78 jieV and 
Ej, H = 82 /ieV. 




FIG. 2: Three- state model of poisoning. Transitions from 
to £ are driven by a nonequilibrium QP source in the leads. 
When AEm < these QPs see the island as a barrier and 
the parity is predominantly even. When AEm > QPs see 
the island as a trap and parity is predominantly odd at low 
temperatures, but QPs can be thermally activated out of the 
island at higher temperatures. 

In Figs. Elb) and (c) we show I sw (n g ) curves for the 
L and H devices measured at T = 30 mK with a ramp 
rate 1 = 1 fiA/s (1 nA/ms). The devices behave quite 
differently: the L device is le periodic while the H device 
is 2e periodic. This behavior has been observed in eigh- 
teen devices made in three different fabrication sessions 
but having very similar values of A^, A^, Eco, and Rn- 
It was also seen in measurements of the zero-bias phase 
diffusion resistance (not shown), where the devices were 
never driven to the voltage state. Since thermal QPs are 
strongly frozen out, the le period in the L device im- 
plies a source of nonequilibrium QPs. Presumably, the 
co-fabricated H device contains a similar QP source, but 
its 2e period seems to indicate that these QPs do not en- 
ter the island. Although this might be naively expected 
from the barrier-like gap profile of the H device, we will 
show below that the actual situation is more complex. 
We begin by constructing a model that explains when 
nonequilibrium QPs in the leads should enter the island 
of an SCPT. 

We consider here three states: the fully paired (no 
QPs) "0" state, the i state with one QP on the island, 
and the 11 £" with one QP in the leads, near a junction. 
We assume QPs cannot be spontaneously created on the 
island, thus poisoning is a two-step process: a QP is cre- 
ated in the leads (transition — > parity remains even) 
and then tunnels onto the island (transition £ — > z, parity 
switches to odd). From the £ state, the QP may also dif- 
fuse away from the junction or recombine with another 
QP (transition £ — > 0), in which case poisoning does not 
occur. For the reduced range n g : — > 1, the energy of 
the state is the ground state energy of the Coulomb 
and Josephson components of the Hamiltonian (see, for 
instance, Ref. [3), E^Tj(n g ,(p), where cp is the super- 
conducting phase difference across the device. The ener- 
gies of the £ and i states are Eg = E^Tj(n g ,cp) + Ai 
and Ei = E^rj(n g ,cp) + A^. The energy change for 
<-» £ transitions is a constant, A^. The energy change 
for parity-switching £ <-> i transitions is 

AEii(ng,(p) = E £ -E % =5E c ^}(n g ^) + 5A. (1) 

Figure [21 illustrates the energy levels of our model for two 
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FIG. 3: (a)-(b) AEu for the L and H devices. The range of 
AEu for all values of <p is indicated by the dark gray region. 
(c)-(d) Predicted crossover m Isw between even and odd par- 
ity. The dotted (dashed) curves are the even (odd) curves 
from Fig. ^a) while the thick lines follow the curve corre- 
sponding to the lowest energy QP state. The island traps 
QPs for rig > n g , C r (light gray regions). 



cases. When AEu < 0, QPs in the leads are prevented 
from reaching the island by a barrier, thus poisoning is 
suppressed. When AEu > 0, QPs will readily tunnel 
onto the island and can become trapped there. This trap- 
ping behavior means that even a weak QP source in the 
leads can result in strong poisoning of an SCPT. 

For | S A | < SEq^j, our model predicts that a single de- 
vice can span both barrier and trap regimes as n g is var- 
ied, regardless of gap profile. Our experiments test this 
prediction quantitatively. During each switching mea- 
surement, rig is fixed but, due to phase diffusion, (p varies 
rapidly as the bias current is ramped |4( . The upper part 
of Figure shows AEu vs. n g with the range of values 
due to variations in <p indicated by a dark gray band. 
The point at which AEu first crosses zero for any value 
of (p defines the critical value n 9yCr above which the sys- 
tem can trap QPs. The lower part of Fig. |3| shows the 
corresponding I sw (n g ) curves predicted by the model as- 
suming the system always occupies the lowest energy QP 
state. For our devices, n^ cr = 0.23 and n^ cr = 0.68 
(with uncertainties of ~20%). 

To compare our data with the predictions above we 
must consider parity fluctuations due to finite temper- 
ature (or external noise) that can cause the system to 
occupy excited states. The ability of our measurement 
to resolve parity fluctuations is determined by the fluctu- 
ation rates Yu(w) for lead-to-island (island-to-lead) QP 
tunneling and the time to ramp between the two switch- 
ing currents at a given n g , r ramp = \(IZ en ) ~ (^)l/^- 
When r rarnp ^> T^jk^ only the lower critical current at 
a given n g is observed since fluctuations into this state 
occur before the ramp can reach the higher critical cur- 
rent. It is important to note that these fluctuations will 



occur even if the occupation probability of the state with 
lower critical current is small. Thus a slow enough ramp 
will only reveal the smaller of if^ en and I°!^ d at a given 
n gi yielding perfect le periodicity as we observe for the 
L device (Fig. [It b)). In Figure 0ta) we shorten r rarnp 
by a factor of 100 and observe that both even and odd 
states are populated, indicating that the parity is indeed 
fluctuating. Since we do not observe both states until 
we ramp at this rate, we estimate YJ^f^ ~ 10 /is — too 
fast to see with our slower ramps. When we show all 
counts for I sw (n g ) equally in Fig. ^b), rather than us- 
ing a grayscale, we see that the odd state is only occu- 
pied when n g > 0.2, in good agreement with the value 
n g,cr ~ 0-23 derived above. 

Parity fluctuations are also apparent in the H device 
when we plot all counts equally (Fig. Etc)). The small 
number of counts below ij^ en indicates occasional fluctu- 
ations to the odd state, occuring for n g > 0.7 as predicted 
Q. This demonstrates that the H device is in fact poi- 
soned, despite the clean 2e curve in Fig.QJc), illustrating 
how a grayscale or average plot of I sw (n g ) can fail to re- 
veal QP poisoning. Comparing Figs.^b) and (c), we see 
that the poisoning is considerably weaker in the H de- 
vice than in the L device. Furthermore, the temperature 
dependence of the two devices is very different: I sw {n g ) 
of the L device changes considerably as T is raised (as 
discussed below) while the H device behavior remains un- 
changed until T > 300 mK when thermal QPs begin to 
populate the island. To explain these aspects of our data 
we return to the model shown in Fig. [21 

Since our model is based on a nonequilibrium QP 
source, we assume the rates r ^o) f° r — > £ (£ — > 0) 
transitions are approximately independent of T. Assum- 
ing £ <-> i transitions are thermally activated, the ratio 
of the corresponding rates is Yu/Yu = exp(— AEu/kT). 
From detailed balance, the steady state ratio between 
the probability of the even (p even = po + pi) and odd 
(Podd = Pi) states is 

a = ^ = ( 1 + ^ ^ = aote~ AE ^ kT , (2) 
Podd \ Yoi J Yu 

where = (1 + Y^/Y^i). Equation predicts the fol- 
lowing for an SCPT with AEu > 0. If there are no QPs 
in the leads, To^ = and the device will always be in the 
even state, as expected. However, for any nonzero To^ the 
device will become trapped in the odd state as T — > 0. 
Furthermore, as T is raised QPs can be thermally acti- 
vated out of the trap, and for T > T* = AEu /kin (aoe) 
the device will be predominantly in the even state (as- 
suming thermal QPs are still negligible). The effect of 
this process on I sw (n g ) is shown schematically in Fig- 
ure Efa). 

With thermal activation, our model explains both L 
and H devices. In the L device I sw {n g ) changes little be- 
low 100 mK, suggesting that the effective temperature of 
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FIG. 4: (a) I sw (n g ) for the L device ramped at / = 100 fiA/s 
and T — 30 mK. The grayscale is saturated at half the max- 
imum count number to emphasize rare events, (b) Same as 
(a) but with all switching events shown equally. The gray re- 
gion marks the predicted range for QP trapping, (c) I sw (n g ) 
for the H device ramped at / = 1 /jlA / s with all counts shown 
equally. (We show data for the slower ramp rate to emphasize 
events that occur below the dominant 2e curve.) 



QPs on the island does not change below this tempera- 
ture. However, above 100 mK, the occupation probabil- 
ity shifts rapidly towards the even state and, at 200 mK, 
the even state is seen more often than the odd state (see 
Fig. [Sfb)). In Fig.[S{c) we display histograms of I sw for 
the L device at n g = 0.7. The histograms show peaks 
corresponding to the even/odd states and their evolu- 
tion with T. To first order we can take the ratio of the 
peak heights, /N!^ x , as a measure of a. The inset 

of Fig. [Sf c) shows that this ratio increases exponentially 
from 100 mK to 200 mK (the peaks become difficult to 
distinguish above 200 mK). A fit using Eqn. yields a 
trap depth AE^f = 62 ± 3 fieV that agrees with the 
value of 72 ± 13 fieV predicted by Equation [I] The fit 
also yields a^f ~ 100, indicating that QPs occupy the £ 
state from which poisoning can occur ~1% of the time. 
We note that behavior consistent with this thermal ac- 
tivation picture was found in at least one other study of 
SCPTs 0. 

Using the prefactor a^f obtained above, we can predict 
the level of poisoning in the H device. For the maximum 
trap depth AEu{n g = 1) = 17 /ieV, Eqn. [21 predicts 
T* = 40 mK. At the minimum effective temperature of 
100 mK reached in the L device, we expect a > 10 at 
rig = 1, i.e., Peven ~ 1- Thus the parity of the H device 
remains mostly even because its gap profile makes AEa 
small and QPs do not become cold enough to remain on 
the island with such a shallow trap. 

In summary, we have demonstrated control over the 
gap profile in SCPT transistors and studied QP poison- 
ing in devices with two types of profiles. We find that 
the behavior of these devices can be quite complex, de- 
pending on gap profile, gate voltage, temperature, and 
measurement timescales. A three- state, nonequilibrium 
QP model correctly predicts the range of gate voltage 
over which QPs can be trapped on the island and the 
temperature dependence of our devices. It also predicts 
that a weak QP source can cause strong poisoning, unless 
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FIG. 5: (a) Thermal activation of QPs out of the island of 
the L device in the trapping regime, (b) I S w(n g ) for the L 
device at 200 mK. (c) Histograms of I sw at n g = 0.7 for the 
L device as a function of temperature. The histograms span 
T — 30 mK to 350 mK in 20 mK increments (from bottom to 
top) and are offset vertically for clarity. Inset: Peak height 
ratio and fit to thermal activation model (see text). 



the trapping regime can be eliminated entirely by making 
a type H device with a larger SA. Finally, our H device 
demonstrates that the traditional method for detecting 
QP poisoning in SCPTs is not completely reliable. Ap- 
plications such as metrology and quantum computing in 
which QP free operation is critical may require better 
tests of QP poisoning, perhaps involving realtime detec- 
tion of individual tunneling events. 

The authors would like to acknowledge useful discus- 
sions with D. Esteve. 
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